Draft version May 13, 2008 

Preprint typeset using MlgX style emulateapj v. 10/09/06 



00 

o 
o 

<N 



Oh- 

o 
is 



(N 

in 
o 

00 

o 



GETTING LUCKY WITH ADAPTIVE OPTICS: 
FAST AO IMAGE SELECTION IN THE VISIBLE WITH A LARGE TELESCOPE 

N.M. LAW^ CD. Mackay^, R.G. Dekany\ M. Ireland^'^, J. P. Lloyd"^, A. M. Moore\ J.G. Robertson^, R Tuthill^, H. Woodruff^ 

Draft version May 13, 2008 

ABSTRACT 

We describe the results from a new instrument which combines Lucky Imaging and Adaptive Optics to 
give the first routine direct diff'raction-Hmited imaging in the visible on a 5m telescope. With fast image 
selection behind the Palomar AO system we obtained Strehl ratios of 5-20% at 700 nm in a typical range of 
seeing conditions, with a median Strehl of approximately 12% when 10% of the input frames are selected. 
At wavelengths around 700 nm the system gave diff'raction-limited 35 milliarc second FWHMs. At 950 nm the 
output Strehl ratio was as high as 36% and at 500 nm the FWHM resolution was as small as 42 milliarcseconds, 
with a low Strehl ratio but resolution improved by factor of ~20 compared to the prevailing seeing. To obtain 
wider fields we also used multiple Lucky-Imaging guide stars in a configuration similar to a ground layer 
adaptive optics system. With eight guide stars but very undersampled data we obtained 300 milliarc second 
resolution across a 30x30 arcsec field of view in i' band. 

Subject headings: instrumentation: high angular resolution — instrumentation: adaptive optics — methods: 
data analysis — techniques: high angular resolution — techniques: image processing 



1. INTRODUCTION 

Adaptive optics has been used successfully for high-angular 
resolution imaging in the near-infrared on a large number of 
5- 10m class telescopes. However, it has not yet demonstrated 
routine diff'raction-limited imaging in the visible on large tele- 
scopes. In this paper we describe the first system capable of 
producing diff'raction-limited moderate- Strehl images in the 
visible on a 5m-cl ass telescope. The system conibines the 
Lucky Imaging (Friel ll978l: [Baldwin et ani2001l: iLaw et al.l 
l2006bh frame selection concept and adaptive optics to pro- 
duce much better performance than provided by either alone. 

Direct Lucky Imaging has demonstrated routine diff'raction- 
limited performance in I-band on 2 . 5m-cla ss telescopes 
(iMackav et al.1 [2QQ4bl: iLiwetaP [2QQ1 l2QQ6allR 120081) . A 
simple, low-cost system has been developed by the Cam- 
bridge Lucky Imaging group and similar systems have now 
been deployed by other groups (eg. Hormuth et al. 2007). 
Lucky Imaging, in contrast to speckle imaging, provides di- 
rect images which can be used irrespective of the complexity 
of the science target. 

We here describe the performance of visible-light Lucky 
Imaging behind an Adaptive Optics (AO) system on a 
large telescope (the experiment is hereafter referred to as 
Lucky-hAO). We used the LAMP (LuckyCam, Adaptive 
optics, aperture Masking and polarimetry at Palomar) in- 
strument, an electron multiplying (EM) CCD-based camera 
custom-built for this experiment to operate in conjunction 
with the Palomar Adaptive Optics system PALM AO. The in- 
strument off'ered three observation modes: Lucky Imaging, 
adaptive-optics assisted aperture masking, and high-contrast 
imaging polarimetry. Results from the latter two modes will 
be described in separate publications. 

During six nights on the Palomar 200" Hale telescope the 
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system routinely produced diffraction-limited resolution im- 
ages with 5-20% Strehl at -700 nm (figure [T]). In this paper 
we characterize the performance of the system under a variety 
of conditions and wavelengths with a view to guide the design 
of more permanent Lucky -h AO instruments. 

In section [2l we describe the observations and instrument 
setup. In section [3] we detail the performance of the system 
and the efliicacy of the combination of Lucky Imaging and AO 
in the low- Strehl regime. We also describe our experiments in 
using multiple guide stars to obtain a larger useful field of 
view. We conclude in section |4l 

2. OBSERVATIONS AND INSTRUMENT SETUP 

LAMP was used in a six night run on the 5m Palomar Hale 
telescope from 2007 July 2 to 2007 July 8. Time equiva- 
lent to two nights was devoted to the Lucky Imaging mode 
described here, during which we observed a variety of both 
test and science targ ets. The Palomar DIMM-MASS unit 
(iThomsen et al.ll2007h was used to record seeing parameters 
during our observations. 

LAMP was placed behind the Palomar Adaptive Optics sys- 
tem PALMAO (Troyetal. 2000; Dekany et al. 1998). The 
system has 241 active actuators with 5yum stroke, and a 16x16 
subaperture Shack-Hartmann wavefront sensor. The system 
operates at up to 2 kHz and typically updated at more than 
200 Hz rates during our observations. In median Palomar 
conditions (V-band seeing =1.1 arcsec) the typical bright-star 
wavefront error is 220nm RMS. We replaced the usual wave- 
front sensor visible/IR dichroic with a 50:50 beamsplitter to 
send visible light to the science focus. 

LAMP was used in place of the AO system's usual near- 
infrared imaging camera. We built a simple reimaging camera 
to provide five field sizes ranging from 31 arcsec (61 milliarc- 
second pixels, undersampled ~4x at 700 nm) to 6.0 arcsec 
(Nyquist sampled at 500 nm). An atmospheric dispersion cor- 
rector was employed to allow broad-band observations. 

The camera detector was based on an electron-multiplying 
528x512 E2V CCD97. The CCD's electron-multiplication 
process allows detection of individual photons in each frame 
produced by the camera at the full quantum efficiency (up to 
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Fig. 1. — Examples of Lucky + AO PSFs acquired with LAMP at 710 nm. Left: The single star HD 160507; a 5% selection from a 600 second run at 20 frames 
per second (linear grayscale). Middle: the same image as the left panel but with a log scale. Right: the binary HD 235089 (Am ~ 4.5, 0.57 arcsec separation) 
observed at 710nm (10% selection from 320 seconds at 50 FPS). 







1 










Fig. 2. — Frames selected from 6000 frames in a typical Lucky Imaging + AO run taken at 50 FPS at 710nm with 14.9 milliarcsecond pixels, targeted at the 
single star HD 192849. Left: one of the high-Strehl frames in the run. Middle: a frame with near-median Strehl. Right: one of the lowest-Strehl frames. 



90%) of the CCD (iMackav et al.ll20Q4al) . Our custom cam- 
era electronics are capable of running at up to 20 frames per 
second (FPS) in full 528x512 pixel frames; subarray readouts 
were used to increase the speed up to 50 FPS for some tar- 
gets. The EMCCD gain was optimized for each target and 
ranged from 1 (no gain) to -10000. The camera produced 14- 
bit data at 7.5 megapixel s per second . The data were recorded 
using custom software (lLawll2007h : a lossless compression 
algorithm reduced the data transfer requirements by an aver- 
age factor of 1 .9, allowing direct recording onto external USB 
hard disks. 

We used the standard Lucky Imaging data reduction 
pipeline without modification; c omple t e desc ripti ons of the 
reduction process can be found in Tlawl (l2007h and Law et alj 
(|2()06b[). The data acquisition software was capable of dis- 
playing a realtime-preview of the Lucky Imaging output but 
the full data reduction was performed by a scripted process 
during daytime operations. Briefly, recorded frames were 
bias-corrected, flat-fielded and cosmic rays were removed. A 
bright star in the field was selected to serve as a Lucky Imag- 
ing guide star (in these observations this was typically the 
adaptive optics guide star). The frames were sorted in or- 
der of Strehl ratio and those that met a specific quality crite- 
rion were selected and aligned to produce a final high-angular- 
resolution image. 

We estimate the Strehl ratio and optimal shift position for 
each image by cross-correlating the instantaneous guide star 
PSF with a diffraction-limited point spread function. The 
height of the peak of the resulting 2D array gives the degree 
of correlation (^ the Strehl ratio) and the position of the peak 
gives the optimal shift position to align the brightest speckles 
of each image. We perform the calculations on an image sub- 
sampled by 4x4 to allow sub-pix el image alignment using the 
Drizzle (iFruchter & Hookll20Q2h algorithm. 
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Fig. 3. — Strehl ratio histograms for three 710nm wavelength observations 
from our dataset with seeings ranging from 0.75 to 1.1 arcsec. Note that for 
clarity the darker histograms obscure parts of the lighter ones. 

We targeted a variety of relatively bright (my = 6-10) 
stars for these observations to ensure that photon noise was 
a small contribution to any observed PSF variability (for our 
faintest observed stars we calculate that photon noise could 
contribute at most | of the observed Strehl ratio variability). 
Frame selection of images containing photon noise can bias 
FWHM measurements of the guide stars ( Law 2007); to avoid 
this we have measured all FWHMs and image profiles using 
the secondary component of close binaries. 

3. PERFORMANCE 

Figure [2] compares three representative frames recorded 
within a single 1 -minute run. It is clear there are large image 
quality variations. Typical Strehl ratio histograms recorded 
by our instrument behind the AO system are shown in figure 
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Fig. 4. — Cuts through a typical Lucky Imaging -i- AO PSF. The companion 
to HD 235089 was observed in 0.86 arcsec seeing at 50 frames per second in a 
10 nm bandpass centered at VlOnm. The y-axis is in data numbers as recorded 
by the camera. The "no correction" (seeing-limited) profile is modeled here 
by a Gaussian with the width measured from the Palomar DIMM during the 
observations and total flux equal to that of our AO-corrected observations. 

O The histograms cover a representative range of seeings and 
frame rates for our dataset; it can be immediately seen that 
the Strehl ratio variations behind the AO system are sufficient 
that large Strehl ratio improvements can be realized by using 
only the frames with the greatest Strehl ratios. 

Each of the Strehl ratio distributions is positively skewed, 
increasing the fraction of high- Strehl outliers compared with 
a Gaussian. The Strehl variability behavior in this low- 
to-moderate Strehl regime contrasts with the negatively- 
skewed distrib utions measured in the high- Strehl regime by 
iGladysz et aP (2006, 2008). Our Strehl ratio distributions 
appear similar to the behavior obse rved for fast-franie-rate 
observations without an AO system (|Tubbs et al.ll2002t iLawl 
120071: [Baldwin etaLll2008h . The Gladysz results were ob- 
tained in K-band on relatively small telescopes, with Strehl 
ratios up to 70% and are thus in a very different ima ge qual- 
ity regime from our data. As noted in [Gladysz et al.l (|2006), 
in the low Strehl regime covering the Lucky -I- AO experiment 
we would expect the observed positively-skewed Strehl ratio 
distributions. 

3.1. Single Star PSFs 

Figure m shows cuts through the output images from a typ- 
ical Lucky Imaging -h AO run. The Airy ring only becomes 
visible in the PSF cuts after frame selection and alignment; 
frame selection leads to improved Strehl ratio and decreased 
FWHM. This is quantified for a typical run in figure [51 Using 
20% of the frames Lucky Imaging reached diff'raction-limited 
resolution for the telescope and provided a >2x improvement 
in FWHM and Strehl ratio compared to using the AO system 
alone. Further improvements in Strehl ratio were achieved 
with more stringent frame selections while the FWHM re- 
mained diffraction-limited. Note that the more stringent selec- 
tions concentrate light from the profile wings visible in figure 
[T] and so appear above the unselected PSF at all radii shown 
in figure [H Figure [6] details the enclosed ffux as a function of 
circular aperture radius for three example Lucky -h AO runs. 

Under the conditions considered here the Strehl ratio in- 
creases rapidly as the quality criterion is made more strin- 
gent. The FWHM, however, is almost diffraction-limited even 
with only shift-and-add alignment and no frame selection (fig- 
ure [5]). Frame alignment reduces the FWHM in size by al- 
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Fig. 5. — Strehl ratio and full-width-at-half-maximum resolution as a func- 
tion of frame selection percentage, for the same run as figure IH The run 
had a duration of 320 seconds during which the minute-timescale seeing was 
0.86 arcsec with a peak-to-peak variability of 0.05 arcsec. AO-only is the 
result of a long-exposure (summing all the recorded frames) behind the AO 
system. 100% "selection" is equivalent to shift-and-add using all the frames. 
Note that the data was slightly spatially undersampled, so the best resolution 
achievable for a diff'raction-limited image is about 35 milliarcseconds. 
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Fig. 6. — The enclosed fraction of total starlight as a function of circular 
aperture radius. Three representative runs are shown (upper lines), all were 
10% Lucky AO selections from runs on HD 192849 and HD 235089 taken 
at 50 FPS in a lOnm bandpass at VlOnm; the seeing was 0.8-1.0 arcsec. The 
dashed line is a seeing-limited profile modeled by a Gaussian with 0.9 arcsec 
FWHM. Note that the structure of all the Lucky -i- AO curves is very similar, 
and the 50% enclosed flux radius is less than 0.2 arcsec. 



most a factor of two, suggesting that a large component of 
the residual wavefront error produced by the AO system is 
in tip/tilt. However, we note that we guide on the brightest 
speckle rather than the centroid of the image, and so in the 
low-Strehl / multiple- speckle regime we obtain much higher 
quality PSFs th an centroiding tip-tilt correction can provide 
(IChristoulll99lh . 

Figure [7] details the performance under a range of seeing 
conditions. In all observations Lucky-hAO significantly im- 
proved on the performance of AO alone. Under median Palo- 
mar Observatory seeing (1.1 arcsec) and selecting 10% of the 
frames Lucky -h AO gave output Strehl ratios of 0.05 to 0.13 
at 700 nm. In 0.7 arcsec seeing the performance increased to 
17% Strehl ratios. In all cases the system gave an improve- 
ment in Strehl of between 2x and 3x compared to the PAL- 
MAO system alone. 
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Fig. 7. — Strehl ratio as a function of seeing. Bottom panel: long exposure 
AO imaging (averaging the frames recorded by camera). Top panel: The 
results of Lucky Imaging of the AO data, selecting 10% of frames. To guide 
the eye, the dashed line in each panel shows the average Strehl ratio achieved 
by the long-exposure AO imaging. To obtain a wide range of seeing the data 
presented here was taken over several hours on two separate nights. Each 
datapoint is a 1.0-4.0 minute dataset taken at frame rates between 20 and 50 
FPS. Data was recorded in 10 nm passbands centered at 670 nm or 710nm; 
the images were sampled at 14.9 mas/pixel. 
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Fig. 8. — Companion detection magnitudes. HD 235089, a 0.57" binary 
observed at 710nm, 50 FPS for 320 seconds is used to calculate the limiting 
contrast ratio at which a companion can be detected with a signal-to-noise 
ratio of 10. The specified binary separation is noted at the right of each line. 
The PSF of the real companion was used to determine the signal within a 
6-pixel aperture. For each radius the RMS noise and its variance was mea- 
sured at 8 apertures distributed around the primary. We conservatively set 
the eff'ective background noise at that radius to be a 2cr excursion above the 
average RMS noise. The bars at the right of the plot give the long-exposure 
performance behind the AO system, using the same data. 

3.2. Companion detection 

The benefits of the results described above are clear for 
crowded field data, where simply increasing the separation 
of stars through high angular resolution is most important. 
In this section we investigate the optimization of the Lucky 
Imaging process for faint, close stellar companion detection. 

As the frame selection criterion is made more stringent the 
light in the PSF becomes more concentrated. This aff'ects 
both the background against which a faint companion must be 
detected and the amplitude of its signal within a diff'raction- 
limited core. However, the smaller eff'ective integration time 
also increases the signal and background noise. The trade- 
off' between these eff'ects is described in figure [8] where we 
show measured companion detection limits for a range of 



close companions vs. the frame selection level for a typical 
observation (based on the binary shown in the right panel of 
figure [T]). 

There is little or no drop in faint companion detection per- 
formance with frame selections as stringent as 10%, at all 
tested radii. At the smallest separations (<0.1 arcsec, not 
shown). Lucky image selection is required to resolve the bi- 
naries. These results demonstrate that, compared to long- 
exposure AO imaging, Lucky-hAO close companion detec- 
tion can be performed with higher angular resolution and with 
minimal SNR performance loss. 

3.3. Performance at other wavelengths 

To investigate the performance of Lucky AO as a function 
of wavelength we observed HD 235089 in a set of bands be- 
tween 500 and 950 nm (figure O. The observations were 
taken during a 45 -minute period during which the DIMM see- 
ing was 0.80+0.15 arcsec. The wavelength range was limited 
in the blue by the lower throughput of the AO system optics, 
and in the red by our CCD detector's sensitivity. As above, we 
use the star's companion (0.6 arcsec separation) to avoid bias- 
ing the PSF shape by the inffuence of photon noise on frame 
selection. 

In all cases the output FWHM resolution was improved by 
a factor of at least two compared to the capabilities of the AO 
system alone. At 500 nm the resolution achieved was as good 
as 43 milliarcseconds, not diff'raction-limited but still a factor 
of 20 improvement over the natural seeing. 

As would be expected, the system gives a decreased output 
Strehl ratio at shorter wavelengths. However, the fractional 
increase in Strehl ratio is actually greater at 500 nm than at 
710 nm, and the frame selection process gives the most obvi- 
ous FWHM improvements at 500 nm. It seems that there is 
increased frame sharpness variance associated with the lower 
quality images at the shorter wavelengths. 

3.4. Frame rates 

Figure[TOlgives Strehl ratio structure functions for three typ- 
ical observations. It is clear from these measurements that at 
our highest speed of 50 FPS the instrument is not completely 
sampling the PSF changes behind the AO system, with most 
of the observed variation occurring on timescales on the level 
of single frames. For stars of the brightness used here photon 
and EMCCD multiplication noise accounts for at most \ of 
the observed single-frame-time variations in Strehl; it appears 
that we really are significantly undersampling the true behind- 
AO Strehl ratio variations. This might be expected given the 
hundreds of frames per second rates typically employed by 
the PALM AO system for turbulence correction. We note that 
future Lucky + AO systems may wish to provide 100+ FPS to 
achieve the highest possible Strehl ratios. 

3.5. Isoplanatic patch 

In the above sections we used the adaptive optics guide star 
as the Lucky Imaging guide star; the greater guide-star sensi- 
tivity of lucky imaging leads the AO guide star to always be 
bright enough for Lucky guiding. If there are multiple guide 
stars in the field, however, we can investigate more complex 
frame quality sensing strategies. In this experiment we use 
multiple guide stars in a 30-arcsecond field in the core of the 
Ml 3 globular cluster. 

To obtain a wide field of view we had to use very spa- 
tially and temporally undersampled data (60 milliarcsecond 
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Fig. 9. — Strehl ratio and F WHM as a function of wavelength and number of frames selected. Each observation was taken at 50 frames per second in a 10 nm 
bandpass filter. As in section [TTI to avoid biasing the FWHM measurement, the companion of HD 235089 was used as the performance measurement star while 
the AO and Lucky Imaging systems guided on the primary. Note that the increased FWHM of the Lucky-Imaging output PSFs at 950 nm compared to 700 nm is 
due to the increased diffraction-limited core size at those wavelengths. Slight undersampling of the original PSFs gives a slightly increased FWHM compared to 
that naively expected for a diffraction-limited PSF on a 5m telescope. 100% selection is equivalent to shift-and-add brightest speckle aUgnment for all frames. 

adaptive optics (GLAO) system. Using eight lucky guide stars 
simultaneously, and selecting the optimal guide star to use for 
each point in the image, the entire 30x30 arcsecond field can 
be processed to give FWHM resolutions better than 300 mas 
in I-band (figure [TT] [right panel]). 
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Fig. 10. — Structure functions, the typical Strehl ratio variance between 
frames as a function of time, for three typical runs taken at 50 FPS at 710 
nm in a 10 nm bandpass. Note the slight downturn at the minimum frame 
spacing of 0.02 seconds, suggesting that we are approaching but not reaching 
the coherence time. 



pixels, 20 FPS). The performance, while much better than the 
prevailing 0.60 arcsec seeing, is therefore greatly decreased 
compared to the results described earlier and what could be 
achieved with a faster, larger array camera system. 

When operating the data reduction system in the usual 
mode of using the AO guide star as the Lucky Imaging refer- 
ence the FWHM approximately doubles over a 15 arcsecond 
radius (figure [TT] left panel). A very diff'erent behavior is ob- 
served when a diff'erent star is used for Lucky Imaging guid- 
ing (figure [TT] middle panel). The performance around the 
Lucky imaging guide star is greatly improved, from 400 mas 
FWHM to 240 mas FWHM, while the performance around 
the AO guide star is reduced. Guiding on the new star caused 
diff'erent frames to be selected as well as diff'erent image shifts 
to be used; the system is not simply correcting tip/tilt aniso- 
planatism. It appears that the AO system occasionally "Luck- 
ily" produced corrections that improved the turbulence expe- 
rienced in other parts of the image. 

The possibility of correcting diff'erent areas of the field with 
diff'erent selections of the same data suggests a method of im- 
proving the isoplanatic patch if multiple guide stars are avail- 
able, in a setup similar to a multiple-guide- star ground layer 



4. DISCUSSION AND CONCLUSIONS 

The combination of Lucky Imaging and Adaptive optics 
produced 35 milliarcsecond resolution at 700 nm on the Palo- 
mar 200" telescope, extending the useful wavelength range of 
the Palomar adaptive optics system into the visible. The per- 
formance detailed in the above sections is unlikely to be the 
best achievable with systems of this type. Our results sug- 
gest that such systems should run very rapidly, as even 50 
FPS observations undersampled the Strehl ratio variations. In 
addition, some of our observations are of lower performance 
than would be expected for the prevailing seeing conditions 
(figure O. This may be due to insufficiently precise non- 
common-path (NCP) error calibrations on some nights. Al- 
though the Palomar AO system has demonstrated NCP cal- 
ibrations as precise as 34 nm wavefront error with the stan- 
dard infrared camera, LAMP typically achieved a Strehl ratio 
of only ~30% on PALMAO's built-in turbulence-free white 
light source in the visible. Future Lucky -h AO instruments will 
likely obtain higher performance as better NCP error calibra- 
tion is achieved. 

Although the observations described here were of relatively 
bright (my = 6-10) stars, we note that the system is quite 
capable of guiding on much fainter stars. Lucky Imaging 
(without an AO system) can guide on stars as faint as 16* 
magnitude (Law et al. 2006b) on 2.5m telescopes and the 
Lucky -h AO technique is therefore capable of using any guide 
star that is bright enough for adaptive-optics operation (as- 
suming broadband imaging and half the light going to the 
AO wavefront sensor). Laser guide star extensions to the 
Lucky -h AO technique will require techniques to mitigate the 
extra background from the laser backscatter, but are possible 
in principle. 

As the first system enabling diff'raction-limited moderate- 
Strehl performance in the visible on 5m-class telescopes, the 
science applications are broad. In particular, crowded field 
visible photometry (for example in extragalactic resolved stel- 
lar population studies) and high-angular-resolution studies of 
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Fig. 11. — Maps of the variation in FWHM over a 30-arcsecond field in the core of the Ml 3 globular cluster, taken with severely undersampled pixels in the 
SDSS i' filter (770 nm central wavelength). The grayscale is set to the same levels in each image for comparison; the variation in FWHM is measured from the 
profiles of 74 stars distributed throughout the field. Left: Lucky Imaging on the AO guide star. Middle: Lucky Imaging on a different guide star. Note the change 
in the variation of the image resolution, becoming elongated between the Lucky Imaging and AO guide stars. Right: an map constructed from the same data 
using 8 different Lucky Imaging guide stars simultaneously and selecting the optimal Lucky guide star for each point in the image. 



nebular emission lines in the visible will benefit greatly from 
the increased resolution. 

The enlarged useful field of view possible with multiple 
guide stars may greatly improve the range and size of fields 
which can be covered. We recommend further investigation of 
the performance of this mode with higher frame rate cameras 
and improved PSF sampling. 

The seeing during our experiments was typical for the Palo- 
mar Observatory site, mostly around one arcsecond. This of- 
fers the intriguing possibility that these resolutions could be 
improved using a larger telescope in better seeing conditions 
- for example, the Keck telescope in 0.5 arcsec seeing may 
be capable of reaching its diff'raction-limit in the visible using 
Lucky + AO techniques. 



The combination of Lucky Imaging and AO off'ers a rela- 
tively simple, low-cost upgrade to current adaptive optics sys- 
tems. As the first direct diff'raction-limited imaging in the vis- 
ible on large telescopes this technique pioneers the science 
which will be performed with future visible- AO systems. 
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